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Abstract. Interpretation of molecular line observations in tenuous circumstellar disks around young G-type stars 
in terms of a disk mass is difficult without a model that describes the chemical structure of these disks. This 
paper now discusses the chemistry in tenuous disks around young solar-type stars based on disk models that take 
into account the presence of a stellar chromosphere. The example of the disk around a 70 Myr old solar-type star 
shows that the dissociating radiation from the chromosphere is stronger than the interstellar ultraviolet radiation 
field up to a distance of ~ 400 AU from the star. Similar to other studies in this research field, it is found that, due 
to photodissociation, the CO-to-H2 ratio is far from the canonical value of 10"^ for molecular clouds. Moreover, 
the dust-to-gas mass ratio, as well as the dust grain size play an important role for the H2 abundance in these 
disks. 
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1. Introduction 

The IRAS and ISO satellites revealed that a large frac- 
tion of nearby pre-main sequence solar-type stars are sur- 
rounded by cool dust. Near- infrared and submillimcter 
observations have shown that the dust is distributed in 
disk-like sometimes ring-like structures around the star 
(see Zuckerman ,2001] for a recent review). However, the 
gas content of these protoplanetary disks is still under de- 
bate. CO observations by Zuckerman et al. H1995|l for a 
sample of young stars in the range 10^-10^ yr lead to the 
conclusion that the gas dissipates very rapidly — within 
a few million years — from the protoplanetary disks. A 
later study by Coulson et al. H1998|l for a sample of 24 
candidate Vega-excess stars shows also too little or no CO 
emission. Only recently, Greaves et al. H2UUU|I found no 
evidence for CO rotational lines in a sample of 14 nearby 
F and G-type stars with known circumstellar dust. Since 
solar-type stars are thought to have little UV radiation, 
the non-detection of CO is generally interpreted as an 
overall lack of gas in the disks around these stars. As a 
consequence planets would have to form either faster in 
the framework of the core-accretion model or they form 
on a very short timescale — a few thousand years — via 
gravitational instabilities. 



What people have neglected so far, is the possible exis- 
tence of a chromosphere in solar- type stars. Ayres II1997|I 
has shown that the ionizing ultraviolet flux from our Sun 
was much stronger in the past and evolved roughly as 
with time. The aim of this paper is to study the influence 
of a chromosphere on the chemistry in late phases of disk 
evolution, that is for disks with Lir/L, < 0.01. For this 
purpose, we choose a template for a 70 Myr old solar-type 
star and compile a typical chromospheric radiation field 
(Sect. 2). The modeling procedure is briefly described in 
Sect. 3 and the outcome of the model calculations is pre- 
sented and discussed in Sect. 4. 



2. The template G-type star 

The initial search for a template star was lead by the fol- 
lowing criteria: 

— The spectral type has to be similar to our Sun, G5V. 

— Lir/L, < 0.01, that is the disk has to be optically thin 
in the UV. 

— Disk parameters have to be known from infrared pho- 
tometry. 

— CO observations have to exist (either detection or up- 
per limit). 
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Table 1. Stellar parameters and disk parameters for the 
template star. The disk parameters in brackets refer to an 
alternative disk model that would result in a very similar 
IR excess. 



parameter template star 



Spectral type 


G5 V 


T'cff 


5570 


logs 


4.5 


stellar radius 


0.9 Rq 


distance 


15.7 pc 


age 


70 10^ yr 


Ri (disk) 


190 AU 




(60 AU) 


Ro (disk) 


470 AU 




4.4 10"^ 


Mdust (disk) 


0.057 M® 




(0.77 Me) 



('1997a' [T997E|l . Zuckerman et al. ifTO^ . Coulson et al. 
(^1998^) and Greaves et al. H2UUU|I . only one candidate star 
was found: HD 123160. Unfortunately, the spectral type 
as well as the nature of the infrared excess of this star was 
recently questioned: Mora et al. I|2()()l|l and Kalas et al. 
H2()()2|l classify this star as a K5 giant and detect a com- 
plex reflection nebula instead of a disk-like structure. We 
conclude that observational data is simply up to now too 
scarce, and we proceed by compiling a template star using 
typical stellar and disk properties. 

2.1. Stellar parameters and disk properties 

The basic stellar parameters and disk parameters assumed 
for the template star are summarized in Tabled We as- 
sume in the following a typical age of 70 Myr for our tem- 
plate. 

Given the fact that observations yield in most cases 
only a lower limit for the gas mass, the disk models are 
calculated for three different gas masses: 0.033, 0.33 and 
3.3 M0. On the other hand, the dust mass is fixed to 
0.057 Mgg. This means that the gas-to-dust mass ratio 
varies in the three models from 0.58 to 58. In addition, 
model calculations are performed with an alternative more 
massive dust disk model, Mdust — 0.77 Mg, which shows 
a very similar IR excess. 

2.2. The UV radiation field 

Following Table 1, a Kurucz ATLAS9 stellar atmosphere 
model (Kurucz 11992(1 with an effective temperature of 
5750 K and a logg of 4.5 is assumed for our template 
star. 

In early stages of star formation, the ultraviolet radi- 
ation field is dominated by active accretion from the cir- 
cumstellar disk/envelope. In later stages, when accretion 



is no longer the dominant process, stellar activity takes 
over. Studies of open star clusters have shown that young 
solar-type stars are more active than the present Sun. 
Ayres l)1997|l estimated correlations between stellar ac- 
tivity and rotation for different tracers like O vi, C ill and 
Ly a based on observed X-ray-rotation and C iv- rotation 
relations. Combining this with observed age-rotation re- 
lations, Ayres concludes that all the photorates scale ap- 
proximately as with time. Photodissociation of CO 
and H2 takes place through discrete bands and is thus af- 
fected mostly by the continuum. Nevertheless, some of the 
dissociating bands overlap with the emission lines of an 
active solar-type star, e.g. with the O vi emission around 
1038 A or the Cm at 977 A. Hence, it is assumed that 
the total flux in the 912 - 1110 A region scales with time 
as ~ t~^. 
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Fig. 1. Flux distribution of the Sun at the solar surface. 
Solid hnc: Kurucz ATLAS9 model with Tcff = 5750 K 
and logg = 4.5. Thin dotted line: LASP/CU rocket ex- 
periment. Thick dotted line: SOLSTICE. 

Data from the LASP/CU rocket experiment and Solar 
Stellar Irradiance Comparison Experiment (SOLSTICE) 
serves as a template for the UV radiation of our present 
Sun. Fig.^shows that the rocket and satellite data agrees 
well with a Kurucz ATL AS9 model of 5750 K and a log g 
of 4.5. 

With an age of 70 Myr for our template star, the flux 
has to be scaled by a factor a = 70 with regard to the Sun. 
Since Ayres l|1997fl states that the linear scaling of the 
flux does not hold for processes involving the continuum 
at wavelength larger than 1500 A, a different scaling law is 
used for the 912 - 1190 A range (rocket data) and the 1190 
- 2000 A range (SOLSTICE data). The first wavelength 
interval is scaled up by a constant factor of 70. To match 
the photospheric flux at 2000 A, the latter is then scaled 
according to the following relation 

log i^i, (template) = log i^^ (SOLSTICE) + 
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Fig. 2. Flux distribution at the stellar surface of the tem- 
plate star. Solid line: Kurucz ATLAS9 model with Toff — 
5750 K and log 5 = 4.5. Thin dotted line: LASP/CU 
rocket experiment data scaled by a factor 70. Thick dotted 
hne: SOLSTICE data scaled according to Eq.lQJ. 

where is given in erg cm~^ s~^ Hz~^ and A is in A. 
The resulting flux distribution is depicted in Fig. [3 Our 
high resolution spectrum overlaps with previously pub- 
lished reconstructed near UV fluxes for the 70 Myr old 
Sun (Dorren & G!uinan ri994|l . The latter does not cover 
the UV region of interest, where CO and H2 photodisso- 
ciate. However, in the overlap region beyond 1250 A, the 
overall agreement is very good. 

Fig. 121 illustrates the strength x of the scaled UV ra- 
diation field with respect to an interstellar Habing field 
(Habing UnSHIl, 1.22210^ cm'^ s'^ x denotes the in- 
tegrated UV flux over the range 912 - 1110 A. It is 
clearly seen that the stellar radiation field dominates up 
to ~ 400 AU in an optically thin disk. 
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Fig. 3. Unattenuated integrated UV flux (912 - 1110 A) 
normalized to the interstellar Habing flux as a function 
of distance from the star: solid line template star plus 
chromosphere, dashed line: Habing field. 



3. The disk model 

The main features of the disk models are briefly summa- 
rized and for further details, we refer to the original paper 
Kamp & Bertoldi (2000). 

We assume thin hydrostatic equilibrium models 

n(r,z) = ni (r/i?i)-2-5 e--V2'»^ _ (2) 

As a simplification, we assume a dimensionless scale- 
height H = hjr ~ 0.15. This does not lead to fully self- 
consistent disk models, but a coupling of the gas tem- 
perature calculation and the disk structure calculation is 
beyond the scope of this paper. The inner and outer ra- 
dius of the disk are set to the respective values given in 
Tabled The power-law exponent of the disk surface den- 
sity is —1.5. The stellar radiation field F^, is described 
above (see Sect. l2.2|l . 

The dust temperature follows from radiative equilib- 
rium assuming large spherical grains of radius a — Z fim 
(Kamp & Bertoldi IMJOjl 

Tdust = 282.5 (L^/Lq)'/' (r/AU)-'/' {a/fim)-'/" , (3) 

with the stellar luminosity in units of the solar luminosity 
Lq . The assumption of this simple form of radiative equi- 
librium is certainly correct for the large dust grains and 
optically thin disks described in this paper. 

The chemical network consists of 47 atomic, ionic, and 
molecular species that are related through 268 gas-phase 
chemical and photoreactions. A number of reactions is 
treated in more detail like H2 and CO photo-dissociation, 
and C ionization. The only surface reactions incorporated 
are H2 formation and freezing out of CO on cold dust-grain 
surfaces. Since we are dealing with large dust particles, we 
reduced the H2 formation rate according to the reduced 
grain surface area. The abundance of CO ice is due to a 
balance between freezing out of gaseous CO and reevap- 
oration of CO ice. The two differences to the Kamp & 
Bertoldi (2000) paper are: (1) the inclusion of cosmic ray 
reactions and (2) a lower temperature of 20 K for freez- 
ing out of CO ice. However, due to the strong stellar UV 
radiation field, the chemistry is mostly driven by stellar 
photons. A modified Newton-Raphson algorithm is used 
to obtain a stationary solution of the entire chemical net- 
work. 

4. Results 

Four main questions are adressed in the following subsec- 
tions: 

— What is the chemical composition of the disk? How 
does it depend on disk mass, what is the role of shield- 
ing? 

— How does the CO-to-H2 ratio change with disk mass? 

— How does the varying dust-to-gas mass ratio affect the 
chemical composition of the disk? 

— What is the effect of a chromosphere as compared to 
the interstellar radiation field? When does a chromo- 
sphere need to be included in disk modeling? 
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Fig. 4. Densities of different species (left scale) and temperature (right scale) in 2 different disk models: (a) 
m0.033_190_CH model at the disk midplane, (b) m3.3.190_CH model at the disk midplane, (c) m0.033_190_CH model 
at one scaleheight, (d) m3.3_190_CH model at one scaleheight. The adopted linestyles are: H:solid, H2: dotted, C: 
dashed, C+: dash-dotted, CO: dash-dot-dotted, COico: dash-dot-dotted (thick), c^: long dashed, T: solid (thick). 



The basic parameters of the disk models computed to an- 
swer these questions are summarized in Table|21 The outer 
radius is fixed to 470 AU in all models. The name of the 
disk model is composed as follows: M^ss-Ri-F^- The type 
of radiation field Fi, is indicated by 'CH' or 'IS', which 
stands for 'chromosphere' or 'interstellar'. The former ra- 
diation field is explained to great detail in Sect. 12.21 The 
interstellar radiation field is taken from Habing (1968) and 
approximated by the following expression 

_ 10-14 / 25 , 25 1 13\ 
" " /I 6 2 3 j 

erg cm"^ s^^ Hz^^ , (4) 

with /5 being the wavenumber in 10^ cm^^. For simplic- 
ity, it is assumed that the interstellar radiation field only 
penetrates the disk from the inner edge. A more realistic 
approach would involve a 2D radiative transfer, which is 
so far not included in the model. 

4.1. The chemical composition of the disk 

We limit the discussion here to some key species, namely 
H, H2, C, C+, CO, and CO ice, that can act as 
possible tracers for the disk mass. The three models 
m0.033_19G_CH, m0.33.190_CH, and m3.3_190_CH cover 
the transition from a completely transparent disk to a 
moderately shielded one. 

Fig. 0] illustrates the chemical structure of two disk 
models: m0.033_190.CH and m3.3_190.CH. Since the disks 



Table 2. Model parameters for the various disk models. 
The columns denote the inner and outer radius in AU, 
the gas mass in Mg, the dust-to-gas mass ratio (5, the dust 
absorption cross section in the UV in cm^ per H-atom and 
the type of UV radiation field. In the last colum, 'CH' 
denotes a pure chromospheric UV radiation field, while 
'IS' stands for an interstellar Habing radiation field with 
Go = l. 



name 


Ri 


Mgas 


5 







m0.033_190_CH 


190 


0.033 


1.73 


3.37(-21) 


CH 


m0.33_190.CH 


190 


0.33 


1.73(-1) 


3.37(-22) 


CH 


m3.3_190_CH 


190 


3.3 


1.73(-2) 


3.37(-23) 


CH 


m0.033.190.IS 


190 


0.033 


1.73 


3.37(-21) 


IS 


m0.33.190.IS 


190 


0.33 


1.73(-1) 


3.37(-22) 


IS 


m3.3.190.IS 


190 


3.3 


1.73(-2) 


3.37(-23) 


IS 


m0.033_60_CH 


60 


0.033 


2.3(1) 


5.38(-20) 


CH 


m3.3_60.CH 


60 


3.3 


2.3(-l) 


5.38(-22) 


CH 


m0.033.60JS 


60 


0.033 


2.3(1) 


5.38(-20) 


IS 


m3.3_60JS 


60 


3.3 


2.3(-l) 


5.38(-22) 


IS 



Note: numbers in parentheses: 3.37(-21) = 3.37 10 



are optically thin in the UV, dust shielding is not impor- 
tant for the photoreactions. H2 can efficiently shield it- 
self against photodissociation within one scaleheight in all 
three models. In the m3.3_190_CH model the molecular 
hydrogen extends even to 1.5 scaleheights. Despite this 
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Fig. 5. The CO-to-H2 ratio e in 4 different disk models: (a) m0.033.190.CH, (b) m0.33.190.CH, (c) m3.3_190.CH, and 
(d) m3.3_190JS. 



high abundance of molecular hydrogen, H2 is not able to 
shield CO against photodissociation. Self-shielding of CO 
sets in for the m0.33_190_CH disk model and becomes effi- 
cient in the highest mass model, m3.3_190_CH. There, the 
CO abundance reaches values of logeco = —5.4. In the 
outermost parts of the disk, r > 390 AU, the temperature 
drops below 20 K, but the low densities ntot < 10^ cm~'^ 
prevent a strong freeze-out of CO on dust grains. 

The chemistry depends very much on the He abun- 
dance: Helium is partially ionized by cosmic rays. The 
ionization degree Heii/Hei is always much smaller than 
1 and depends strongly on the recombination processes 
and hence on the densities in the disk model. Since we 
use here the solar He abundance of logenc = 0.1, the de- 
structive reaction He+ -I- CO keeps the C+ abundance up 
throughout the disk in all models. 

4.2. The C0-to-H2 ratio 

Fig. |5l illustrates how the The C0-to-H2 ratio changes 
throughout the disk for various disk models. With in- 
creasing disk mass, the CO abundance in the disk rises 
and hence the CO-to-H2 ratio increases. In general, the 
difference between a model with chromospheric radiation 
field and one with interstellar radiation field amounts at 
maximum to A CO-to-H2 ratio ^ 0.5 dex. 

For an interpretation of observations, which give rather 
an integrated picture of the CO in the disk, it is better 



Table 3. [CO/H2] ratio for a sample of disk models, 
name [C0/H2] 



m0.33_190_CH 


3.58 10" 


7 


m3.3_190.CH 


1.12 10" 


5 


m0.33.190.IS 


1.89 10" 


6 


m3.3.190.IS 


6.22 10" 


5 


m3.3_60.CH 


7.53 10" 


5 


m3.3.60JS 


3.65 10" 


4 



to use the integrated [CO/H2] ratio of the various disk 
models (see Table 



[CO/H2] = 



Je{CO)dV 
I<ii2)dV 



(5) 



Given the tenuous character of these disks, an IS UV 
radiation field would penetrate the disk from all directions 
and therefore lead to a much higher [H2/CO] ratio even for 
the 3.3 model. Thus, we draw the tentative conclusion 
that the [CO/H2] conversion factor in these disks can be 
orders of magnitude smaller than the canonical value for 
molecular clouds, [C0/H2]~ 10"'*. These disks resemble 
in their nature generally the photon dominated regions 
(PDR) on the surfaces of molecular clouds, very high UV 
irradiation and densities ranging from 10'^ to a few times 
10^ cm-3. 
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Fig. 6. Densities of different species (left scale) and temperature (right scale) in 2 different disk models: (a) m3.3_60_CH 
model at the disk midplane, (b) m3.3_60JS model at the disk midplane, (c) m3.3_60_CH model at one scaleheight, 
(d) m3.3_60JS model at one scaleheight. The adopted linestyles are: H:solid, H2: dotted, C: dashed, C"*": dash-dotted, 
CO: dash-dot-dotted, COico: dash-dot-dotted (thick), c~: long dashed, T: solid (thick). 



4.3. The role of dust 

As already stated above, dust shielding plays no role for 
the photochemistry in the disk. On the other hand, the 
dust grain surfaces are the location where molecular hy- 
drogen is formed. Hence, the H2 formation rate 



R„. = 3 X 10-17 ( cm^s-i (6) 

^ VlOOK J 6 X 10-22cm2 ^ ' 

is higher in models with a higher dust-to-gas mass ratio, 
e.g. m0.033_190_CH as compared to m3.3_190_CH. cr„^, 
denotes here the effective UV attenuation cross section 
per H nucleus. Assuming a low dust-to-gas mass ratio of 
1.7310^2 for the 0.033 model prevents formation of 
H2 throughout the whole disk. 

The alternative model (Sylvester & Skinner ri996|l with 
an inner radius of 60 AU and a dust mass of 0.77 has 
an even higher dust-to-gas mass ratio (m3.3_60_CH: 0.23 
compared to m3.3_190_CH: 1.7310^2) and therefore H2 
extends up to two scaleheights and CO up to one scale- 
height in the disk. Additionally, the freezing out of CO is 
slightly enhanced, but still not efficient enough to obtain 
a high CO ice abundance. But even in these disk models 
with a high dust mass, dust absorption is not a dominant 
shielding process. 



4.4. The role of the chromosphere 

The outer parts of the disk beyond 400 AU are clearly 
dominated by the interstellar radiation field. Nevertheless, 
for the inner gas disk the chromospheric UV radiation is 
the driving mechanism in chemistry. If one neglects the 
presence of a chromosphere, the resulting H and CO abun- 
dances differ by more than a factor 10, especially in regions 
where most of the disk mass is located (see Fig. EJ. 

For somewhat younger stars, the chromosphere will be 
even more important, because the stronger stellar activity 
will give rise to an even higher UV photon flux. At a cer- 
tain point back in time, the disk has been optically thick 
and actively accreting; therefore UV photons and X-rays 
from the accretion process have probably determined the 
chemistry and temperature balance in the upper layers of 
the protoplanetary disk. At which stellar age the transi- 
tion to optically thin disks occurs as well as the point in 
time when accretion ceases and stellar activity takes over 
has still to be determined. 

5. Conclusion 

In this paper, we investigate the chemical structure of ten- 
uous circumstellar disks around young solar-type stars us- 
ing the disk models developed by Kamp & Bertoldi H2000|l . 
Previous observational studies have neglected the possible 
influence of a chromosphere for solar-type stars. Especially 
in an earlier evolutionary stage (WTTS), the UV radia- 
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tion from the chromosphere is much stronger than that of 
the present Sun and may determine the chemistry in the 

circumstellar surrounding. 

The main results of this work are: 

— We present a recipe to calculate the chromospheric UV 
radiation field of young solar-type stars as a function 

of age. 

— The CO-to-H2 ratio in low mass disks around solar- 
type stars is smaller than the canonic value of 10~^ 
for molecular clouds and resembles more the values 
found for the PDR surface layers of molecular clouds. 

— The dust-to-gas mass ratio has a strong influence on 
the molecular hydrogen abundance in disks. Low mass 
disks, in which hydrogen would otherwise be atomic, 
can become molecular in the presence of a high dust- 
to-gas mass ratio. Note that since these disk models 
are optically thin in the UV, this is not due to dust 
shielding, but due to an enhanced H2-formation rate. 

— The chromosphere of a 70 Myr old star dissociates the 
CO in the inner disk (r < 400 AU) of a solar-type star. 
Due to self-shielding effects, H2 is much less affected 
than CO. In the outer disk, the stellar UV radiation 
field is smaller than the typical interstellar UV field. 

This work presents a first step towards the modeling of 
protoplanetary disks around young solar-type stars. The 
model is for example limited by the fac;t that we did not in- 
clude an interstellar UV radiation field in a self-consistent 
way (2D radiative transfer). Moreover, at least a simple 
line radiative transfer has to be included in order to (a) be 
able to calculate the gas temperature in these disks self- 
consistently and (b) to extend these disk models to higher 
masses. 
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